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FIELD  DISTRIBUTIONS  IN  THIN  WATER  FILMS - 
A  COMPARISON  OF  THEORETICAL  AND  EMPIRICAL  TRENDS 

Dawn  L.  Scovell,  Tim  D.  Pinkerton,  Valentin  K.  Medvedev,  and  Eric  M.  Stuve 

University  of  Washington 
Department  of  Chemical  Engineering 
Box  351750 

Seattle,  WA  98195-1750 


ABSTRACT 

Field  ionization  ot  thin  water  layers  adsorbed  onto  a  platinum  field 
emitter  was  investigated  by  numerical  simulation  and  analysis  of 
experimental  data.  The  numerical  simulation,  which  includes  a  field- 
dependent  relative  permittivity,  was  developed  to  predict  the  field 
distribution  around  a  water-covered  field  emitter  tip.  The  model  predicts 
that  the  dominant  field  occurs  at  the  water-vacuum  interface  in  thin 
layers.  In  the  experiments,  water  adlayers  were  grown  under  field-free 
conditions  by  exposure  of  a  cryogenically  cooled  emitter  tip  to  water 
vapor  in  ultrahigh  vacuum.  Field  ionization  was  probed  by  ramped  field 
desorption  (RFD)  in  which  desorption  of  ionic  species  (hydrated 
protons)  is  measured  while  increasing  the  applied  electric  field  linearly 
in  time.  The  onset  field  of  ionization  decreased  from  0.3  to  0.2  V/A  as 
temperature  increased  from  130  to  150  K.  An  activation  barrier  of  0.7 
eV  (16  kcal/mol)  for  ionization  of  water  to  produce  hydrated  protons  and 
hydroxide  ions  was  estimated  from  the  temperature  dependence  of  the 
onset  field.  The  experimental  trends  agree  with  the  predicted  trends  for 
thin  water  layers. 


INTRODUCTION 

High  surface  electric  fields  drive  many  important  processes,  such  as  » 
electrochemistry,  corrosion,  and  field  emission.  These  surface  fields,  which  are  of  the 
order  of  magnitude  of  1  V/A,  are  strong  enough  to  make  and  break  chemical  bonds. 

The  behavior  of  water  in  these  fields  is  important  because  water  is  the  primary 
component  in  electrochemical  and  corrosion  processes  as  well  as  a  major  contaminant 
in  the  vacuum  surrounding  the  field  emitter  arrays  in  flat  panel  displays.  It  is  usually 
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assumed  that  water  amplifies  the  fields  at  electrode  surfaces  (1),  but  little  is  actually 
known  about  how  water  affects  the  electric  field  distribution. 

Field  emitters  lend  themselves  to  the  study  of  the  dielectric  properties  of  water 
because  they  allow  independent  control  of  the  electric  field  for  fields  as  high  as  5  V/A. 
Multilayers  of  water  can  be  vapor  deposited  by  exposing  a  cryogenically  cooled  emitter 
tip  to  water  vapor.  The  water  deposited  under  these  conditions  is  either  amorphous  ice, 
crystalline  ice,  or  a  mixture  thereof.  Under  the  appropriate  conditions  amorphous 
water-ice  can  undergo  large  translational  movement,  however,  crystalline  layers  are 
stationary  (2-4). 

In '  this  paper  we  report  numerical  and  experimental  results  for  the  field 
ionization  of  water.  ,  Water  was  adsorbed  at  temperatures  of  130  to  150  K  onto  a 
platinum  field  emitter  tip  under  field-free  conditions  in  ultrahigh  vacuum.  This 
temperature  range  was  selected  because  the  deposited  water  layer  is  crystalline  and  will 
be  immobile  during  the  course  of  the  experiment.  Ionization  was  examined  by 
isothermal  ramped  field  desorption  (RFD)  performed  as  a  function  of  temperature  and 
water  layer  thickness.  We  interpret  the  results  in  terms  of  the  thickness  dependence  of 
the  field  at  ionization  onset  and  compare  the  experimental  results  to  those  of  the 
numerical  model  (5). 


NUMERICAL  MODEL  RESULTS 

The  applied  field  required  to  ionize  a  water  layer  increases  with  increasing 
thickness.  This  is  easily  explained  by  the  dielectric  screening  of  water;  thicker  water 
layers  are  more  effective  at  screening  the  field  at  the  water-vacuum  interface.  In  the 
numerical  simulation  the  water/tip  system  was  modeled  in  spherical  geometry.  Key 
assumptions  were  that  the  tip  is  a  perfect  conductor  and  that  no  ions  are  present  in  the 
water  layer  (5).  The  potential  distribution  was  found  by  numerically  solving  Poisson’s 
equation  for  the  water  layer  and  Laplace’s  equation  for  the  vacuum.  The  relative 
permittivity  of  the  water  was  allowed  to  vary  between  low  and  high  field  limits  of  80 
and  2,  respectively,  with  the  transition  occurring  between  0.01  and  1  V/A.  The 
potentials  in  the  water  layer  and  vacuum  were  calculated  as  a  function  of  tip  potential 
and  water  thickness. 

The  model  predicts  that  the  field  is  always  concentrated  at  the  tip-water  and 
water-vacuum  interfaces.  Fig.  1  shows  the  dependence  of  constant  electric  field  lines  '  > 
on  water  thickness  and  applied  electric  field.  The  predominant  field  is  predicted  to  be  at 
the  water-vacuum  interface  in  thin  water  layers.  As  the  water  thickness  increases,  the 
field  at  the  vacuum  interface  decreases  and  the  field  at  the  tip  surface  increases.  When 
the  water  becomes  sufficiently  thick,  the  location  of  the  predominant  field  shifts  to  the 
tip-water  interface.  The  left-hand  curve  in  Fig.  1  represents  the  locus  of  points  where 
the  field  at  the  water-vacuum  interface  Fv  is  0.4  V/A  while  the  right-hand  curve 
represents  the  same  for  the  field  at  the  tip-water  interface  Ft.  Ionization  begins  when 
either  Fv  or  Ft  becomes  sufficiently  large.  As  illustrated  in  Fig.  1,  ionization  is  expected 
to  begin  at  the  vacuum  interface  in  thin  layers  and  at  the  tip  surface  in  thick  layers. 
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Positive  ions  formed  at  the  vacuum  interface  will  be  ejected  from  the  water  layer 
when  they  are  created,  as  shown  in  Fig.  2a.  This  implies  that  the  onset  of  ionization  is 
not  influenced  by  the  negative  ions  that  are  left  in  the  adlayer.  Positive  ions  formed  at 
the  tip  surface,  on  the  other  hand,  must  migrate  through  the  adlayer  and  then  desorb 
from  the  surface  before  they  can  be  detected,  as  shown  in  Fig.  2b.  Thus,  ions  will  be 
present  in  the  water  layer  before  they  are  detected  in  our  measurements.  A  key 
assumption  in  the  numerical  model  is  that  no  ions  are  present  in  the  adlayer.  Clearly 
this  is  not  true  when  ionization  occurs  at  the  tip  surface  and,  as  a  result,  the  numerical 
model  is  not  expected  to  accurately  predict  the  onset  of  ionization  in  thick  layers. 

EQUIPMENT  AND  PROCEDURE 

The  experimental  apparatus  and  procedure  have  been  described  previously  (6) 
and  will  therefore  be  described  only  briefly.  The  vacuum  chamber  was  pumped  by  a 
turbo-molecular  pump  and  titanium  getter  pump  to  obtain  a  base  pressure  of  10  10  torr. 
A  field  emitter  tip  was  used  to  generate  electric  fields  as  large  as  4  V/A.  Ions  produced 
by  the  tip  were  projected  onto  a  pair  of  chevron  microchannel  plates  (MCP)  and 
directed  onto  a  phosphor  screen  to  form  a  greatly  magnified,  projection  image  of  the  tip 
surface.  This  image  shows  the  detailed  atomic  arrangements  of  the  hemispherical  tip 
surface  (7).  Field  ionization  images  were  recorded  with  a  CCD  camera  and  VCR. 

The  tip  was  spot-welded  to  a  0.25  mm  diameter  heating  loop  mounted  to  the 
bottom  of  a  cryogenically  cooled  downtube.  The  temperature  was  measured  by  a 
chromel/alumel  thermocouple  spot-welded  to  the  heating  loop.  The  desired  temperature 
was  maintained  within  ±0.3  K  by  a  PID  temperature  controller  connected  to  a  d.c. 
power  supply  for  resistive  heating.  A  high  voltage  feedthrough  at  the  bottom  of  the 
downtube  provided  high  voltage  isolation  for  the  heating  loop/tip  assembly. 

The  field  emitter  tip  was  electrochemically  etched  from  0.013  mm  diameter 
platinum  wire  at  2-3  V  d.c.  in  a  molten  mixture  of  sodium  chloride  and  sodium  nitrate. 
The  tip  radius  was  determined  from  the  best  imaging  field  of  3.75  V/A  for  neon  (8). 
The  relationship  between  the  applied  potential  V,  and  the  applied  field  Fapp  is 


Fapp  ~  Vf  !  P  rr  > 


[13  ' 


where  (3  is  the  shape  factor  and  r,  is  the  tip  radius.  The  shape  factor  is  5  for  a  typical 
emitter  tip  (9).  The  tip  radius  in  these  experiments  ranged  from  335  to  374  A.  Tip  radii 
determined  from  the  best  imaging  field  were  within  20%  of  those  estimated  by  counting 
the  number  of  lattice  steps  between  the  (001)  and  (113)  planes  in  the  field  ionization 
images.  The  relative  error  in  the  electric  fields  reported  in  these  experiments  is  less 
than  5%.  In  this  paper  the  term  “applied  field”  is  used  to  describe  the  field  that  would 
occur  for  bare  tip  surface  in  vacuum  for  a  given  applied  tip  potential. 
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A  three-step  procedure  was  used  to  ensure  that  the  tip  was  free  from 
contaminants  at  the  beginning  of  each  experimental  run.  First,  the  tip  was  heated  to 
500  K  to  desorb  any  water  from  the  tip  and  heating  loop  assembly.  Second,  any 
remaining  contaminants  were  field  desorbed  from  the  tip  by  imaging  in  neon  and 
hydrogen.  Finally,  a  small  amount  of  platinum  was  field  desorbed  in  neon  to  ensure 
that  the  surface  was  atomically  clean. 

Prior  to  dosing,  water  was  degassed  with  both  liquid  nitrogen  and  dry  ice  in 
several  freeze-thaw  cycles.  Water  was  introduced  into  the  background  of  the  vacuum 
chamber  through  a  variable  leak  valve  to  produce  a  uniform  layer  of  ice  on  the  tip  under 
field-free  conditions.  The  ti^  was  held  at  the  desired  temperature  during  dosing.  The 
tip  was  exposed  to  5.0  x  10  torr  of  water  for  2  to  60  minutes.  Pressures  of  less  than 
5.0  x  10  torr  were  typically  achieved  within  2  to  10  minutes  after  dosing.  The  ice 
layer  thickness  was  estimated  from  the  molecular  flux  during  dosing  and  adjusted  for 
any  desorption  during  chamber  pump  down.  The  sticking  coefficient  was  assumed  to 
be  unity,  and  the  ice  density  was  varied  linearly  from  0.8  g/cm3  at  80  K  to  0.93  g/cm3  at 
130  K  (10).  Water  thicknesses  of  50  to  3600  A  were  obtained. 

Ramped  field  desorption  (RFD)  spectra  were  obtained  after  removing  water 
from  the  background  gas  (1 1).  During  the  RFD  runs,  the  potential  applied  to  the  tip  was 
increased  at  a  constant  rate.  When  the  field  becomes  sufficiently  large,  ions  are  created, 
desorb  from  the  water  layer,  radially  accelerate  away  from  the  tip  and  are  then  detected 
when  they  arrived  at  the  MCP/phosphor  screen  assembly.  The  signal  from  the  phosphor 
screen  was  monitored  by  two  rate  meters,  each  with  a  time  constant  of  0.1  s.  The  first 
rate  meter  was  always  set  to  a  sensitivity  of  300  counts  per  second  (CPS)  full  scale, 
while  the  second  was  adjusted  to  capture  the  entire  peak.  The  first  rate  meter  was  used 
to  determine  ionization  onset. 

EXPERIMENTAL  RESULTS 

The  effects  of  temperature  and  ice  thickness  on  the  ionization  of  water  were 
Studied.  The  ionization  onset  field  Fappo  defined  experimentally  as  the  applied  field  at 
which  ionization  reaches  10  counts/s,  was  measured  as  a  function  of  the  dimensionless 
ice  thickness  at  temperatures  between  130  and  150  K.  The  dimensionless  thickness  x  is 
the  ratio  of  the  thickness  of  the  ice  layer  to  the  tip  radius.  The  ionization  threshold  of 
10  counts/s  is  lower  than  used  previously  (6)  because  the  background  noise  levels  were 
lower  in  these  experiments.  The  applied  field  was  ramped  at  6  mVA  V1  for  all  runs. 

Fig.  3  shows  that  the  dependence  of  the  onset  field  at  133  K  for  thin  layers  of  water  is 
linear  with  thickness. 

The  slope  dF app.o! dx  is  strongly  dependent  upon  temperature.  Fig.  4  shows  the 
temperature  dependence  of  dFapp_0/dx  for  thin  water  adlayers.  The  open  circles  are  the 
experimental  data.  The  slope,  d Fapp,0  Idx ,  ranges  from  0.2-0.3  V/ A.  The  slopes  and 
error  bars  were  determined  by  linear  regression.  The  solid  line  in  Fig.  4  is  the  slope 
predicted  by  the  numerical  model. 
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Extrapolation  of  the  data  in  Fig.  3  to  zero  thickness  gives  an  estimate  of  the 
actual  Field  needed  for  ionization,  F0.  Extrapolation  is  necessary  to  correct  for 
screening  of  the  Field  by  water  layers  of  Finite  thickness.  The  onset  Field  F0  is  a  true, 
rather  than  applied,  Field  because  the  influence  of  water  layer  thickness  has  been 
removed.  The  extrapolated  onset  Fields  are  shown  in  Fig.  5.  The  trend  in  the 
extrapolated  onset  Field  is  linear  between  130  and  150  K. 

DISCUSSION 

The  numerical  model  predicts  that  ionization  occurs  at  the  vacuum  interface  in 
thin  water  layers.  Fig.  1  shows  that  in  the  limit  of  zero  thickness,  the  applied  field  is  the 
same  as  the  Field  at  which  water  ionizes  at  the  vacuum  interface  Fv.  Thus,  the  vacuum 
field  is  the  same  as  the  onset  Field  Fa.  From  Fig.  1,  FV  is  0.31  V/A  at  133  K.  Values  for 
Fv  as  a  function  of  temperature  are  shown  in  Fig.  5. 

The  solid  line  in  Fig.  3  is  the  predicted  trend  for  Fv  =  0.31  V/A.  The  predicted 
trend  agrees  well  with  the  trend  of  the  experimental  data.  The  best  way  to  evaluate  the 
numerical  model  is  to  compare  the  predicted  and  measured  values  of  dFapp,0/dx  because 
the  value  of  the  vacuum  field  (i.e.,  the  y-intercept  of  Fig.  3)  was  determined  from  the 
experimental  data.  The  temperature  dependence  of  the  predicted  slope  is  the  solid  line 
in  Fig.  4.  The  temperature  dependence  of  the  predicted  slope  is  within  the  error  of  the 
experimentally  obtained  slopes. 

The  temperature  dependence  of  the  onset  field  of  ionization  provides 
information  about  the  energetics  ot  the  ionization  event.  As  discussed  in  a  previous 
publication  (6)  the  rate  of  thermally  excited  ionization  may  be  written  as: 


rion 


/ivv\)exp 


V\F„+Ha-{tJ>)F„ 

kT 


[2] 


where  nw  is  the  areal  density  of  water,  v  the  pre-exponential  factor,  V,  the  field 
dependence  of  the  potential  of  curve  crossing  between  the  neutral  and  ionic  potential 
energy  curves,  F0  the  onset  field,  Ha  the  field-free  activation  barrier,  A P  the  difference 
between  the  dipole  moment  of  the  final  and  initial  states,  k  Boltzmann’s  constant,  and  T 
temperature.  This  expression  assumes  that  the  neutral  and  ionic  potential  energy  curves 
are  straight  in  the  region  of  crossing. 

As  a  constant  detection  rate  was  used  to  determine  the  onset  of  ionization,  an 
instrumental  constant  Ic  can  be  estimated  from  Eqn.  2: 
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I  c  =-ln 


rion,o 

nwv 


kT 


[3] 


where  rlono  is  derived  from  the  specified  detection  rate  of  10  counts/s  and  the 
hemispherical  area  of  the  tip.  The  quantity  in  square  brackets  in  Eqn.  2  is  assumed  to  be 
independent  of  field.  The  instrumental  constant  is  45  for  the  specified  detection  rate 
and  a  tip  radius  of  350  A  (6).  Solving  for  F0  gives 


1 

V,  -(& P) 


(H„-IckT). 


[4] 


This  expression  predicts  a  linear  relationship  between  the  onset  field  for  ionization  and 
the  temperature.  As  shown  in  Fig.  5  the  RFD  data  are  linear.  The  linear  relationship  is 
in  excellent  agreement  with  results  reported  by  Stintz  and  Panitz  (12). 

Extrapolation  of  the  data  to  zero  field  provides  an  estimate  of  the  field  free 
activation  barrier.  As  shown  in  Fig.  5,  the  extrapolated  temperature  for  crystalline  ice  is 
216  K,  which  corresponds  to  a  field-free  activation  barrier  for  ionization  of  0.7  ±  0.1 
eV.  This  compares  well  with  results  from  d.c.  conductivity  experiments  for  water 
between  platinum  electrodes  which  show  that  0.74  eV  is  required  to  form  a  pair  of  ions 
from  two  solvated  water  molecules  (13).  The  error  in  the  field-free  activation  energy  is 
large  because  the  data  were  taken  over  a  relatively  narrow  temperature  range  and  have 
large  error  bars.  It  may  be  possible  to  obtain  a  better  estimate  of  the  field-free 
activation  energy  for  the  crystalline  phase  by  extending  the  temperature  range.  This  can 
be  done  by  depositing  crystalline  ice  at  temperatures  above  145  K  (14,  15)  and 
performing  the  RFD  experiment  at  lower  temperatures. 

The  slope  in  Fig.  5  can  be  used  to  estimate  the  dipole  moment  of  the  layer.  If  V; 
is  negligible  (6)  then  the  slope  of -0.003614  VA^K'1  is  equivalent  to  a  dipole  moment 
of  -5.1  Debye.  If  the  final  state  hydronium  ion  is  assumed  to  have  a  small  dipole 
moment  (P/  =  0),  then  this  result  gives  a  value  of  5.1  Debye  for  the  initial  state  water 
molecule.  The  dipole  moment  of  liquid  water  has  been  calculated  to  be  2.4  -  3.0  Debye 
and  that  of  hexagonal  ice  is  predicted  to  be  3.1  Debye  (16).  A  likely  explanation  for  the 
differences  between  the  experimental  and  calculated  values  is  that  V;  is  not  zero. 
Unfortunately,  the  relative  importance  of  the  curve  crossing  (Vy)  and  the  dipolar 
interaction  effects  (AP)  cannot  be  measured  with  this  method,  and  further  research  is 
necessary  to  devise  a  means  for  measuring  these  quantities  independently. 
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CONCLUSIONS 

The  ionization  of  water  layers  adsorbed  under  field-free  conditions  depends 
upon  the  thickness  of  the  water.  Our  numerical  model  predicts  that  thin  water  layers 
ionize  at  the  vacuum  interface.  The  trends  measured  in  our  ramped  field  desorption 
experiments  are  consistent  with  the  model  predictions  for  thin  layers.  The  experimental 
data  were  used  to  estimate  a  field-free  activation  barrier  of  0.7  eV,  which  agrees  well 
with  the  value  of  0.74  eV  obtained  by  others.  The  experimental  data  was  also  used  to 
estimate  a  dipole  moment  of  the  layer  of  5.1  Debye. 
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Fig.  1.  Lines  of  constant  vacuum  field  Fv  and  tip  field  F,  of  0.4  V/A  as  a  function  of 
dimensionless  thickness  x.  Dimensionless  thickness  is  the  ratio  of  the  water 
thickness  to  the  tip  radius.  The  dashed  line  is  the  thickness  at  which  the 
predominant  field  changes  from  the  water-vacuum  to  the  tip-water  interface.  In  thin 
layers  ionization  is  expected  to  occur  at  the  vacuum  interface  and  in  thick  layers  it 
is  expected  to  occur  at  the  tip  surface. 


Fig.  2.  Schematic  of  the  ionization  process  at  (a)  the  water-vacuum  interface  and  (b)  at 
the  tip  surface.  The  heavy  lines  represent  the  interface  where  ionization  occurs. 
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Fig.  3.  Applied  field  at  ionization  onset  as  a  function  of  deposited  ice  thickness  x  at 
133  K.  The  solid  line  is  the  numerical  prediction  for  a  vacuum  field  Fv  of  0.31 
V/A.  The  circles  are  the  experimental  data. 


Fig.  4.  Temperature  dependence  of  the  slope  dFapPt0  /d.t  for  thin  coverages.  The  open 
circles  are  the  measured  data  and  the  solid  line  is  the  numerical  prediction. 
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T/K 

Fig.  5.  The  temperature  dependence  of  the  onset  field  F0  required  to  ionize  thin  layers  of 
water.  This  field  was  determined  by  the  y-intercept  in  plots  of  applied  field  Fappo 
vs.  thickness  .t  (e.g.,  Fig.  3).  Note  that  in  the  numerical  prediction  F0  is  equivalent 
to  the  field  at  the  vacuum  interface  Fv. 


